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NYO-9193 
Abstract 

An  attempt  is  described  to  confine  a  plasma  composed  of 
electrons  and  Ions  by  static,  one-component ,  electric  and 
magnetic  fields.   All  quantities  depend  on  one  space  coordinate 
only,  and  the  particle  acceleration  vector  has  two  components. 
The  contained  plasma  is  assumed  to  be  at  rest  with  no  electric 
or  magnetic  fields  present.   For  the  major  portion  of  the 
report,  the  electron  and  ion  temperatures  are  assumed  equal 
and  no  trapped  particles  are  permitted.   For  both  species  of 
particles  the  Vlasov  equation  is  used. 

A  charge-neutral  approximation  is  discussed  together  with 
the  exact  solution,  and  numerical  results  are  presented.   It  is 
shown  on  one  hand  that  under  the  given  assumptions  the  boundary 
layer  cannot  be  charge -neutral  throughout,  and  on  the  other 
hand  that  even  the  exact  solution  does  not  correspond  to  a 
desirable   confinement  configuration,  as  it  contains  a  runaway 
ion  jet  and  large  electric  fields. 

Finally,  two  simple,  albeit  impractical,  cases  are  pre- 
sented, one  with  unequal  ion  and  electron  temperatures,  and  the 
other  with  trapped  particles,  in  which  a  containment  situation 
is  indeed  obtained. 
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The  Boundary  Layer  In  the  Confinement  of  a 
One -Dimensional  Plasma 

I .   Introduction 

In  this  report  we  describe  attempts  to  obtain  a  simple 
configuration  in  which  a  plasma  composed  of  both  ions  and 
electrons  is  contained  by  static  electric  and  magnetic  fields. 
We  consider  the  simplest  non-trivial  situation  In  which  the  one- 
component  fields  depend  on  one  space  coordinate  only  and  the 
particle  acceleration  vector  has  only  two  non-vanishing  compo- 
nents.  We  assume  the  contained  plasma  to  be  at  rest  and  in 
thermal  equilibrium  with  equal  ion  and  electron  temperatures. 
We  also  assume  that  no  fields  are  present  in  the  contained 
plasma.   We  require  that  there  are  no  trapped  particles  in  the 
boundary  layer  between  the  confined  plasma  and  the  vacuum;  that 
is,  we  demand  that  all  particle  orbits  start  and  end  in  the  con- 
fined plasma.   We  assume  that  the  particles  interact  with  each 
other  only  through  the  average  field  produced  by  all  the  parti- 
cles; that  is,  we  employ  the  Vlasov  equation  for  both  species. 

The  next  section  contains  the  explicit  formulation  of 
the  problem  and  a  detailed  description  of  the  trapped 
particles.   We  discuss  the  "charge-neutral"  approximation  and 
then  give  certain  expansions  useful  in  the  numerical  treatment 
of  this  approximation  and  the  exact  solution.   The  results  of 
a  numerical  computation  of  the  charge -neutral  curve  are  in 


1.  For  earlier  work  directly  related  to  this  problem  but  with- 
out electric  field,  see  H.  Grad,  Phys.  Fluids  4,  1366  (1961) 
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the  third  section.   We  conclude  that  a  boundary  layer  without 
trapped  particles  cannot  te  essentially  charge  neutral  through- 
out the  layer.   In  the  fourth  section  we  present  the  numerical 
solution  to  the  full  problem.   We  find  that  the  only  solution 
Is  physically  unacceptable  as  It  contains  a  run-away  Ion  jet 
and  large  electric  fields.   We  are  thus  forced  to  consider 
more  complicated  containment  configurations  and  In  the  last 
section  we  present  two  simple  schemes,  one  with  different 
electron  and  Ion  temperatures,  and  the  other  with  trapped 
particles.   We  offer  these  examples  only  to  show  that  some  kind 
of  confinement  Is  Indeed  possible;  we  do  not  propose  them  as 
physically  Interesting  --  or  attainable  --  states. 


II.   Formulation  of  the  Problem 

Consistent  with  the  conditions  of  the  previous  section, 
we  assume  that  the  electric  and  magnetic  fields  vary  with  the 
X  coordinate  only,  that  the  electric  field  has  an  x  component 
only,  that  the  magnetic  field  has  a  z  component,  and  that  the 
non-trlvlal  velocity  components  are  u  In  the  x  direction  and 
V  In  the  y  direction.   The  Vlasov  equation  for  the  distribu- 
tion function  of  each  species  of  particles  (+  Ions,  -  electrons) 
becomes 


Sf+         e^  (  af^     af+" 


(1)       u  ^  (x,u,v)  +  —   (E  +  vB)  ^  -  uB  ^ 


=  0 


-  R  _ 


and  the  fields  satisfy  Maxwell's  equations  which,  in  rationalized 
MKS  units  are 


:2a) 


K 


11^  =  e(n_^(x)  -  n_(x)) 


(2b) 


-^i^=e(j.(-)  -J-W)  ' 


where 


(3a) 


and 


:3b) 


n+(x)  = 


f ,  (x,-u,v)dudv  , 


J+(x) 


'f^(x,u,v)dudv 


In  order  to  solve  the  Vlasov  equation  and  to  intro- 
duce the  concept  of  trapped  particles,  we  define  the  character- 
istic coordinates  of  (l)  by 


(^a: 


52LM  =  u(s 
ds  ^ 


(4b) 


^^-  I   (E(x(s))    +   v(s)B(x(s))) 


;4c) 


dv(s)   ^ 
ds 


I  u(s)B(x(s))    , 


so  that 


(5) 


^  f(x(s),    u(s),    v(s))    =   0 
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Equation  (5)  states  the  well  known  result  that  f(x,u,v)  Is 
constant  on  the  trajectories  of  the  particles  In  the  given 
fields.   There  are  two  sets  of  equations  (4)  and  (5)>  one  set 
for  the  ions  and  one  for  the  electrons;  however,  we  tempor- 
arily suppress  the  distinction.   We  know  that  the  solution 
of  the  equations  (4)  has  the  exact  constants  of  the  motion 

(6a)  ^  ^  2  ^^^     "*"  ^  )  +  e^ 

(6b)  P  =  mv  +  eA  , 

where  we  have  introduced  the  usual  scalar  and  vector  potentials 
by 

(7a)  ^li^  =  -E(x) 

dx 

(7b)  ^  =  +B(x)  . 

In  a  certain  sense  we  might  say  that  a  class  of  solutions  to 
equation  (5)  is  any  function  f(e,P). 

In  the  problem  we  are  considering  we  assume  the  distribu- 
tion function  to  be  Maxwellian  as  x  — >  -00,  or 

n  m        ,    /  2,   2x 

„/  N  O  ,  1  (U   +   V    )  N 

f(x,u,v)  =  2B^  exp(-  2  m  -^-^ ^)  . 

X  — >    -00 

As  we  suppose  that  $(x)  and  A(x)  tend  to  zero  as  x  tends  to 
minus  infinity,  we  see  that 
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n  m 
(8)  f(x^u,v)  =  2^j^  exp(-  ^) 


One  might  Infer  from  the  above  remarks  that  (8)  holds  every- 
where In  the  plasma,  but  this  need  not  be  the  case.   Certain 
trajectories,  defined  by  {^)  ,    need  not  reach  x  =  -oo;  thus  the 
"boundary  condition"  (8)  will  not  determine  f(x,u,v)  on  those 
orbits.   We  cail  such  particles  trapped  particles,  and  through- 
out the  main  portion  of  this  work  we  assume  that  no  trapped 
particles  are  present.   Prom  (6)  we  see  that 


2  ?  2 

m  u   =  2m(£  -  e$)  -  (P  -  eA) 


so  that  for  u   to  be  real. 


(9)  2m(e  -  e$)  -  (P  -  eA)^  >  0 


Given  a  particle  at  a  point   x  with  energy   e   and  canonical 
momentum  In  the  y  direction   P,  the  above  Inequality  must  hold 
at  all  points  to  the  left  of  Its  position  (i.e.,  for  all 
x'  <  x)  If  the  particle  Is  to  be  able  to  reach  (or  to  come 
from)  the  point  at  -oo,  or.  In  other  words,  for  the  particle 
not  to  be  trapped.   For  such  untrapped  particles  we  assume  (8) 
holds;  otherwise  we  take  f(x,u,v)  =  0. 

We  may  now  easily  eliminate  equation  (l)  entirely  by 
expressing  n^(x)  and  J^(x)  In  terms  of  A(x)  and  |(x).   After  a 
change  of  variables  from  u   and  v   to   e   and  P  we  obtain 
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from  the  basic  definitions  (3)j  the  equations 


(10a)     n(x) 


n 


o 


TTkT 


dedP  e 


-e/kT 


J  J 

Dfx 


(10b)      J(x)  = 


n 


(  f 


o 


niTrkT 


y2m(e-e^(x))  -  (P-eA(x))' 
dedP  e"^'^^'^  (P-eA(x)) 


j^|4)y2m(£-e$(x))  -  (P-eA(x))' 


where  the  domain  of  integration  D(x)  is  defined  by  the  restric- 
tion that  it  contains  all  values  of  e,   0  <  e  <  oo  and  all 
values  of  P,   -oo  <  p  <  oo,  such  that 

(11)       2m(e-e$(x')  -  (P-eA (x ' ) )^  >  0 

for  all  x'  <_  x.   In  order  to  proceed  further  we  must  recognize 
the  difference  between  the  ions  and  electrons.   For  later 
reference  we  non-dimensionallze  (lO)  and  (11 )  fur  ions  and 
electrons  separately. 

We  Introduce  non-dimensional  potentials  by 


(12a) 


U{^) 


e+  f(x) 


kT. 


+ 


(12b) 


a^(x)  = 


e,A(x) 


/2m^kT 


+ 


On  performing  the  e  integration  firsts  and  with  the  definitions 
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c 


(13)      F(x)  =  e^ 


and 


e  '^  dy 


X 


(14) 


e+(p.x)  -   sup{;5+(x')  +  (p-a+(x'))  \   , 
x'<x 


we  obtain 


(15a)     n+(x)  = 


2n  (^ 

o 


I  \ ^  -e+(p,x) 

dp  P  ye+(p,x)-j25  (x)-(p-a+(x))^  le  " 


(15b)     J+(x)  = 


2"o  \¥^ 


vim 


-+ 


dp(p-a, (x))e 


-e, (p,x) 


?(/£+(?, x)-j2^+(x)-(p-a+(x))' 

the  non-dimensional  variable  of  Integration  p  being  defined 
as   P//2krrL|.T_^"  respectively  for  the  Ion  and  electron  case.   The 
system  of  equations  Is  now  (2)  and  (15). 

In  order  to  complete  the  non-dlmenslonallzatlon  of  the 
full  system,  we  define 


(l6a) 


.1/^ 


(l6b)      T  =  T   (=  T   by  assumption) 


"+ 


(l6c)     a(x)  =  -ya.    (x) 


-  a, (x) 


-  10 


(I6d)  j6{x)   =  A  (x)   =   -t>_{^ 


+ 


(l6e) 


1   \\  ^^^     , 


and  the   equations   become 


(17a; 


p+{i)    ==  I  i   dp  G+{vA)e      - 


(17b)     j+(|)  =  I 


r 


dp  G+{vA)e      -  (p  +  Y   a(| 


(I8a) 


^a(^)  =  ^  j_(U  -  TJ+(^) 
d^         ' 


(18b)       ^^%p  ^i^(4)  -  P_(e)  -  P^(^), 
(m_m_|_)  ^   d^ 


where 


(19a)     e+(p,U  = 


sup  s+jZ5(|')  +  (p  +  Y"^  a(^'))^> 


i^<i 


and 


TiT77;2 


(19b)     G+(p,U  -  F  /e+(pT)77(^)^  (P  +  Y'^  a(4 


The  dimensionless  constant  in  (l8b)  is  just  ^    (Vc)  ,  where 
V   Is  the  geometrical  mean  of  the  ion  and  electron  thermal 
speeds  and   c   is  the  speed  of  light. 
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Before  solving  the  full  system  of  differential  equations 
we  consider  an  approximation  of  both  mathematical  and  physi- 
cal meaning.   Since  v/c  should  be  much  less  than  one,  it  is 
reasonable  to  assume  that  the  solution  of  (l8b)  is  given 
approximately  by  the  solution  to 

(20)       0  =  p_^(4)  -  p_(|)  . 

Equation  (20)  states  that  the  plasma  is  charge-neutral,  and 
physically  we  expect  that  the  plasma  should  try  to  remain  in 
an  approximately  charge-neutral  state.   Mathematically  we 
expect  that  the  solution  to  (l8)  could  be  developed  in  an 
asymptotic  expansion  in  powers  of  KM-kT//m_m   and  that  the 
leading  term  would  be  given  by  the  solution  to  (20).   Since 
there  is  no  explicit  ^  dependence  in  (20),  we  may  reinterpret 
the  equation  as  a  definition  of  an  implicit  functional  rela- 
tion between  f)   and  a.   With  this  functional  relation  known 
we  might  solve  the  one  remaining  differential  equation  (l8a). 
We  shall  obtain  by  numerical  means  the  functional  relation 
between  ^   and  a   for  various  values  of  y.   We  do  not  proceed 
to  the  solution  of  the  remaining  differential  equation  as  we 
would  not  obtain  any  further  information  of  interest. 

In  order  to  check  and  to  start  the  numerical  work  we 
have  obtained  an  expansion  of  the  solution  to  (20),  the  charge- 
neutral  curve,  for  a  and  ^   small.   After  relatively  straight- 
forward manipulations  we  derive: 
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(21) 


i^ 


-i  a  +  ma 
o 


5/2  +  na^A  +  o(a2)  , 


where  i   Is  the  unique  non-negative  zero  of 


(22a) 


f  (i)  =  /2l 


/x  dx  e 


-(x 


i  ,2    ,-° 


J 

o 


2Y^v^ 


/^  dx  e 


-(x-#)' 


ana 


(22b) 


and 


m 


o 


TT 


f'(i  )  <  ° 

^  o 


(22c) 


n 


9' 


In  formulas  (21),  (22),  we  assume  y  <  1  (as  it  is  physically). 
One  may  also  verify  as  easily  that  there  is  no  other  root  of 
(20)  of  the  form  fi   -  ia"  +  ...  ,  n  >  0,  starting  from 
j?5  =  a  =  0.   In  the  next  section  we  give  the  full  numerical 
solution  to  (20). 

In  order  to  begin  the  solution  to  the  full  differential 

equations  (l8a)  and  (l8b),  we  need  starting  values  of  ^ ,   -^  , 

Sa 
a,  -^   ;    moreover,  we  must  show  that  there  is  a  solution 

leaving  the  singular  point  j?5  =  a  =  0  of  the  equations  (l8). 

If  we  look  for  a  solution  which,  for  small  ^   and   a,   has  the 

expansion 


(23) 


fi   =  -ia  +  . 


(i  >  0)  , 
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then  we   obtain   from    (l8a)    and    (l8b)    the   expansions 


,2 


(24a)  ^  ^(|)    ^  a(i)  /irm  + 


d^' 


(24b)  -£ 


where 


(25) 


/" 


a(i)    =  ^{(2t)^/^ 


o 


/x    (x-^)dx  e 


-(--2^ 


2 


-<!' 


3/2 


XI 


dx 


/^   (x-4?)    e 


-(x-4?)^ 


J 


and  f(^)  Is  defined  by  (22a).   From  (24a)  we  infer  that 


(26) 


a(^)  =  (^)'(4-^,)V 


^  >  I 


-  0  '    ^  <  ^o  • 


As  we  may  choose  ^   arbitrarily,  we  set  I   =0.   If  we  Insert 
this  expansion  in  (24b)  we  obtain  the  relation 


(27)      f(i)  - 


-  f  ^a(i) 


KM-kT 
/m_m 


For  Kli.kT//m_m   sufficiently  small  we  infer  the  existence  of 
a  root  of  (27)  near  the  value  i   at  which  f(£    )  =  0.   Thus  we 
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have  shown  the  existence  of  a  solution  to  the  full  differential 
equations  which  leaves  the  singular  point  jZ5  =  a  =  0.   The  form 
of  this  solution  is  given  by  the  relations  (23),  (26),  and  (27). 
We  do  not  assert  that  this  solution  is  unique;  however,  there 
is  no  other  solution  of  the  type  (23),  and  with  some  effort  one 
can  prove  that  there  is  no  solution  of  the  form  ;Z5  =  "t  ia  +  .  .  .  , 
n  >  0.   Thus  we  have  chosen  to  ignore  some  solutions  in  which 
electric  effects  are  larger  than  magnetic  ones  and  badly  behaved 
solutions,  e.g.,  those  with  oscillations  of  period  tending  to 
zero  as   a   and  p   tend  to  zero.   We  note  further  that  the  solu- 
tion obtained  is  in  fact  near  the  solution  given  by  the  charge- 
neutral  curve.   In  the  fourth  section  we  describe  the  numerical 
solution  of  (l8) . 


III.   The  Charge-Neutral  Curve 

Before  describing  the  results  of  the  computation  of  the 
solution  of  (20),  we  explain  briefly  the  techniques  used.   The 
function  P(x)  appearing  in  the  definition  of  the  charge  and 
current  densities  satisfies  the  differential  equation 

(28a)      F'  (x)  =:  -1  +  2xF(x) 

and  has  the  integral  representation 

00 

r 


;28b)      F(x)  ^  \ 


e   du 
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We  computed  F(x)  for  one  large  value  of  x  from  (28b)  and 
then  Integrated  (28a)  Inward  to  set  up  a  table  of  values  of 
F(x).   Since,  for  x  >  0,  we  are  integrating  in  a  stable 
direction  we  expected  no  errors  to  propagate  and  we  found 
that  the  error  in  the  table  of  values  was  no  larger  than 

_o 

a   few  times  10   .   Besides  tabulating  F(x),  we  also  stored 
the  first  three  derivatives  of  F(x),  so  that  we  could  use 
Taylor's  formula  for  interpolation.   Thus  the  total  error  in 

_o 

the  computation  of  F(x)  was  no  larger  than  a   few   times  10 
The  only  real  complication  in  the  computation  of  the  solution 
of  (20)  comes  from  the  definition  of   e+(p,x)   by  (I9a).   We 
stored  tables  of  values  of  the  variable  of  integration  p   in 
(l7a)  and  the  associated  values  of   £+(p,x).   After  solving 
for  a(x)  for  a  given  )6{yi) ,   we  corrected  the  values  of  £_).(p,x) 
according  to  (l9a),  and  then  we  increased  ^{x)    and  searched 
for  a  solution  of  p   -  p_  =  0.   In  the  searching  procedure 
the  first  guess  was  a  linear  extrapolation  of  the  previous 
results  or  the  value  given  by  (21).   If  this  value  was  not 
sufficiently  close  to  the  root,  then  we  tried  another  and 
finally  interpolated  among  the  various  values  until  a  suffi- 
ciently precise  root  was  found.   We  discovered  that  unless 
p   -  p_  was  quite  small  (not  much  larger  than  10~')  cumulative 
errors  could  destroy  the  solution.   Further,  we  found  that  in 
critical  sections  of  the  curve  a  relatively  small  increment  in 
^{x)    each  step  was  also  needed.   Except  for  these  difficulties 
the  numerical  computation  was  quite  direct.   We  might  add, 
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however,  that  the  computing  time  for  the  charge-neutral  curve 
was  almost  as  long  as  the  time  to  obtain  the  solution  to  the 
full  differential  equations. 

We  computed  at  least  part  of  the  charge-neutral  curve 
for  three  mass  ratios,   !Ti_/m  =  -,  o^/  ,  y^^  ,  -=r-Q    •   '^^®  compu- 
tations were  easiest  and  most  reliable  for  the  smallest  mass 
ratio.   In  Figure  1  we  plot  ^   and  n  (=  n_)  as  a  function  of 
a   for  the  mass  ratio  ttt^  -■   The  curves  for  the  other  mass 
ratios  were  qualitatively  the  same  as  In  Figure  1,  but  while 
the  curves  turned  back,  we  could  not  be  sure  of  the  details  of 
the  turn-around  since  the  finite  net  of  points  necessary  for 
the  Integrals  (l7a)  limited  the  ultimate  accuracy  In  this  very 
sensitive  region  of  the  curve.   In  another  connection  In  the 

next  section  we  give  the  corresponding  curve  for  a  mass  ratio 

„  _1_ 
°^   100  • 

We  considered  earlier  the  possibility  of  expanding  the 
solution  of  the  full  differential  equations  in  an  asymptotic 
series  in  the  parameter   K[J-kT//m_m  ,  and  the  first  term  of 
this  series  was  to  come  from  the  charge-neutral  solution.   We 
infer  from  Figure  1  the  impossibility  of  such  an  expansion 
throughout  the  full  boundary  layer,  for 


(^5)    -  =  -i-iii-iB 


and  we  see  that  bj6/^a   is  Infinite  at  the  turn-around  point. 
Thus  we  are  forced  to  turn  to  a  numerical  solution  of  the  full 
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differential  equations  in  order  to  discover  what  happens  to 
the  solution  in  this  region. 


IV.   The  Solution  to  the  Full  Differential  Equations 

In  order  to  solve  the  system  of  differential  equations 
(18),  we  used  the  methods  of  the  previous  section  to  compute 
p(x)  and  j(x)  and  then  the  modified  Runge-Kutta  method  to 
advance  the  solution  one  step.   Then  we  adjusted  the  quanti- 
ties  s_,_(p,x)  to  include  the  new  values  and  proceeded  with 
the  integration.   We  checked  the  results  primarily  by  com- 
paring with  the  charge -neutral  solution  and  also  by  recomputing 
then  with  different  Initial  conditions  and  different  meshes. 
With  the  meshes  chosen  none  of  the  quantities  changed  more 
than  a  few  percent  in  any  interval  (except  at  the  very  start 
of  the  solution  where  its  behavior  is  given  by  (26)).   We 
shall  examine  several  characteristics  of  the  solutions  obtained. 

We  checked  first  that  the  solution  we  obtained  did  in 
fact  follow  the  charge-neutral  curve  initially.   In  Figure  2 
we  plot  j25  as  a  function  of  a  for  m_/m^  =  j^   for  several 
values  of  the  ratio  v/c;  we  also  Include  the  charge  neutral 
curve  for  comparison.   In  Figure  2  we  see  exactly  the  behavior 
we  would  expect  if  the  charge  neutral  curve  is  to  be  the  leading 
term  of  an  asymptotic  expansion  of  the  full  solution:   as  the 
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expansion  parameter  becomes  smaller  and  smaller  the  solution 
approaches  the  charge- neutral  curve. 

In  figures  3?  ^,    and  5  we  present  the  results  of  the 
numerical  computations  for  various  values  of  the  parameters 
involved.   In  Figure  J> ,      m  /m  =  Too'  ^/.'^  l/s)  =  -p^y  ;  in 
Figure  4 ,   m_/m  =  Too  ^  ^A*^  ^^^  "^  Tn  '  ^^^   ^^   Figure  5  >  we 
chose   m_/m  =  -.  o^/  and  v/c  /sj  =  -ftt.   We  have  chosen  the 
peculiar  variable 

(30)       E(x)  =  E(x)    ^^^^ 

2)/m_m 

for  then  the  electromagnetic  part  of  the  pressure  tensor  assumes 
the  simple  form 


(31)      pelectromagn.^  1  ^^2  _  ^2^ 


Thus  we  see  that  after  the  solution  leaves  the  charge- neutral 
curve  the  electrons  turn  back  very  rapidly  and  only  the  ions 
remain.   The  resulting  net  positive  charge  causes  the  electric 
field  to  increase  as   x   increases,  and  this  in  turn  acceler- 
ates the  ions  forward.   Since  the  magnetic  field  still  turns 
back  particles  the  ion  density  drops.   However,  the  ion  pres- 
sure must  continue  to  rise  to  balance  the  electric  field 
pressure.   In  other  words,  we  find  evidence  of  a  type  of  run- 
away phenomenon,  with  a  jet  of  highly  energetic  ions  leaving 
the  plasmao   This  solution  obviously  is  not  acceptable  physically 
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as  it  requires  electric  fields  whose  pressure  is  much  larger 
than  the  magnetic  pressure.   Thus  the  conclusion  that  we 
offer  is  that  there  is  no  steady-state  confinement  with  equal 
ion  andelectron  temperature  without  trapped  particles  in  the 
boundary  layer.   In  the  next  section  we  show  that  by  relaxing 
either  of  two  of  the  conditions  imposed  we  can  find  a  confine- 
ment situation. 


V.   Other  Solutions. 

If  we  relax  the  condition  that  both  species  of  particles 

must  be  at  the  same  temperature,  then  we  may  easily  obtain  a 

solution  with  j25  =  E  =  0.   If  we  consider  a  situation  in  which 


(31)       T_/T_^  =  m_^/m_ 


then  we  see  from  (12)  and  (15)  that,  since  a  (x)  =  -a_(x) 
the  assumption 


(32)      J25  -  E  -  0 


is  consistent  with  the  necessary  condition,  n  (x)  =  n_(x). 
Hence  the  equations  reduce  to  the  magnetic  ones  alone  and  the 
solution  becomes  the  one  found  by  H.  Grad  ,  after  a  simple 

T"!   See  footnote,  page  4. 
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rescallng  of  parameters  appropriate  to  a  system  with  mass 


(33a )     —  = 1-  — 


and  temperature 


(33b)     T  =  T^  +  T 


We  may  find  another  type  of  confinement  for  a  plasma 
with  equal  Ion  and  electron  temperatures  by  the  addition  of 
trapped  particles.  .Again  we  look  for  a  solution  with 
j25  =:  E  =  0.   Hence  we  must  dem.and  that  P,  =  P_  everywhere. 
We  note  that  since  there  are  no  electric  fields  zero  velocity 
particles  are  stationary.   Hence  we  can  add  an  arbitrary  density 
of  trapped  particles  of  zero  velocity  at  any  point  in  the  plasma. 
In  particular,  we  may  consider  the  introduction  of  zero-velocity 
electrons,  in  order  to, match  the  net  positive  charge  of  the 
free  particles.   Since  the  trapped  electrons  introduced  are  at 
rest,  they  do  not  contribute  to  the  current;  hence  equation 
(l8a)  still  holds.   The  condition  of  charge  neutrality  obviously 
determines  the  density  of  the  trapped  electrons,  and 

trapped/  \  _   free/  n    „f^ee 


(3^)      P_  ^^   (x)  -  p_^   \yi)    -   p_   (x)  . 


1 


Details  of  such  a  solution  are  shown  in  Figure  6  for  —  -   Yon' 
We  see  that  the  electrons  turn  back  rapidly  and  that  the  ions 
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turn  back  much  more  slowly.   The  characteristic  size  of  such 
a  layer  is  thus  the  ion  larmor  radius.   We  do  not  claim 
physical  relevance  for  these  solutions;  we  offer  them  only  to 
show  that  the  problem  of  confinement  is  really  a  question  of 
picking  the  correct  solution  out  of  many  possible  ones. 
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Figure  I.    Charge  neutrality  curve  for  mass   ratio  =  1/1836. 
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Figure  2.   Dependance  of  the  solutiori  on  v/c  for  the   case 
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Figure  3a.     Solution  for  m_/m^  =  1/100;    v/c/2  =  1/50. 
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Figure  3b.    Solution  for  m_/m^=  l/lOO;  v/c/2  =  1/50. 
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Figure  4a.    Solution  for  m_/m+ =  1/ 100;    v/ct^  =  I/IO. 
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Figure  4b.       Solution  for  m_/m^  =  l/lOO;     v/c/2  =  l/IO- 
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Figure  5a.    Solution  for  m./m^=  1/1830,'   v/c/2  =  1/50. 
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Figure  5b.     Solution  for  m_/m^=   1/1836;   v/c/2  =    1/50. 


1.8 


-  30   - 


KjU> 


in 


m 


ro 


-    51    - 


OCT  2 

1962 

DATE  DUE 

MAR  4    U 

AUG   s    -^ 

GA YLORO 

=  «    -,'  ZD     •*   ^     S    A 

rTYU  c.    X 

NYO-Q193 

Pasklevici 

The  boundary  layer  in  the  con- 
finement of  a  one-dimensional 
plasma. 


N.  Y.  U.  Courant  Institute  of 
Mathematical  Sciences 

4  Washington  Place 
New  York  3,  N.  Y. 


